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ABSTRACT 

We present a study of strong cool core, highly-luminous (most with 
Lx ^ 10^^ erg s~^ ), clusters of galaxies in which the mean jet power of the 
central active galactic nucleus must be very high yet no central point X-ray source 
is detected. Using the unique spatial resolution of Chandra^ a sample of 13 clusters 
is analysed, including A1835, A2204, and one of the most massive cool core clusters, 
RXCJ1504. 1-0248. All of the central galaxies host a radio source, indicating an active 
nucleus, and no obvious X-ray point source. For all clusters in the sample, the nucleus 
has an X-ray bolometric luminosity below 2 per cent of that of the entire cluster. 
Most have a nucleus 2 — 10 keV X-ray luminosity less than about 10^^ erg s~^ . 
We investigate how these clusters can have such strong X-ray luminosities, short 
radiative cooling-times of the inner intracluster gas requiring strong energy feedback 
to counterbalance that cooling, and yet have such radiatively-inefficient cores. If the 
central black holes have masses ~ lO^M© then the power exceeds one per cent of 
the Eddington luminosity, and they are expected to be radiatively-efficient. Only 
if they are ultramassive (Mbh > lO^^M©), would their behaviour resemble that of 
lower mass accreting black holes. Our study focuses on deriving the nucleus X-ray 
properties of the clusters as defined in the above question, while briefly addressing 
possible solutions. 

Key words: Galaxies: clusters: general - X-rays: galaxies: clusters - cooling flows - 
galaxies: jets 



1 INTRODUCTION 

Many clusters of galaxies have steeply rising X-ray surface 
■ brightness profiles, drops in temperature by a factor of 2 to 
\ 3 nea r the centre (|Allen et al.ll200ll : IPe Grandi k Molendil 
I2OO2I ) and short radiative cooling times. These are known 
as cool core clusters. With peaked X-ray surface bright- 
ness, cool core clusters cool rapidly in their centres. To 
maintain pressure support, this results in an inflow of mat- 
ter on the scale of hundreds of solar masses per year, 
also known as a cooling-flow. However, both Chandra and 
XMM — Newton show much less cooling gas below a fac- 
tor of 2 to 3 of the outer temperature than is predicted from 
the c o oling flow mode l (Bohringer et al. 2001; Tamura et alj 
i] 12001. ,2003 : Peterson fc Fabian ,20061 : 
McNamara k Nulsenl l2 



Peterson et"all|2001 . ,2003 : ,Peterson fc Fabian ,20061 : 
20071 ). 

Additional heating is needed to offset the cooling. For 
cool core clusters, the central active galactic nucleus (AGN) 
is likely to be the energy source heating the gas by interact- 
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in^ with the ICM ([Birzan et aLllJOOi : iRaffertv et"alll2006l : 
iDunn k Fab"ia3l2006l . 120081 )^ 

Energy from the AGN can be injected into the ICM 
either radiatively or mechanically. Extrem e examples of the 
first are quasars and bright Seyferts fsee iHopkins k, Elvid 
[2OIOI, for an example on how radiation can couple via a wind 
with the surrounding medium), whereas the second can be 
seen through outflows and jets in radio galaxies. Cygnus A is 
an ideal example hosting a very powerful nucleus. Chandra 
images reveal that this cluster has an obvious point-like core, 
and a spectral analysis shows that the absorption-corrected 
core X-ray luminosity accounts for more than 50 per cent 
of the cluster X-ray luminosity. Radio images of this source 
revea l large scale jets interacting with the ICM (| Young et al.l 
I2OO2I ). 

Many authors have shown in recent years that AGN 
can energetically offset the cooling in galaxy clusters 
by inflating large cavities correlated with radio lobes, 
and propagat ing energy through sou n d /pressure waves 
(iFabian etaP l2003l: feirzan et all 1 20041: [Ruszkowski et"ar 



20041: iForman et al.1 l2005l: [Punn fc FabianI [JOOel: 



Fabian et al.1 I2OO6I : IRaffertv et al.1 [20061 : ISiiacki k Springell 
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Figure 1. Chandra images from left to right of R1504, A1835, R1347 and Cygnus A, with one pixel (e. g. 0.492^^) binning. The colour 
scale is linear expect for Cygnus A which is in a log scale. Black shows regions with zero counts. North is up and East is to the left. The 
top panels are in the 0.5 — 7 keV energy range, and the bottom panels are in the 5 — 7 keV energy range. The X-ray images of Cygnus A 
clearly reveal its point-like core, whereas those of R1504 and A1835 do not show any obvious evidence of one. Also shown are the images 
of R1347 which begin to show a hint of point-like core. 




2QQ6I: iM cNamara k NulsenI l2007l : ISanders fc FabianI l2007l : 



Dunn fc Fabian 2008). However, not all cool core clusters 
have an obvious point-like core in their X-ray emission like 
that of Cygnus A. 



iMerloni Sz Heinj (|2007l ) studied a sample of 15 sub- 
Eddington AGN for which there was an estimate of the ki- 
netic jet power available from the literature, and a clear 
nuclear X-ray source reported (although 2 sources only had 
upper limits of their nuclear source). From there, these au- 
thors were able to find a correlation between kinetic power 
and nuclear X-ray luminosity. Our study presents a sample 
of strong cool core, highly- luminous ( Lx ^ 10^^ erg s~^ ) 
clusters of galaxies, for which their is no evidence of a nu- 
clear X-ray source in the Chandra images. We select a sam- 
ple with known, luminous, cool cores from the Chandra 
archives. It includes some of the most extreme objects such 
as RXCJ1504.1-0248 fsee lOgrean et al.ll201Ql ). All of the 
clusters in our sample require very energetic feedback from 
their AGN to offset cooling, and all have a radio source, in- 
dicating that the AGN should be switched on, yet none have 
an obvious point source in the X-ray band. Unless the central 
black hole is ultramassive (Mbh ^ lO^M©), our objects are 
operating at high enough Eddington rates that they should 
be radiatively efficient. This makes t he lack of any evident X - 
ray nucleus quite puzzling (see also iMcNamara et al.ll20ld V 



In Section 2, we present the sample and reduction tech- 
niques. In Section 3, the results are shown and in Section 4, 
they are discussed. We adopt Ho — 70 km s~^ Mpc~^ with 
Qm = 0.3, Qa = 0.7 throughout this paper. 



2 SAMPLE SELECTION AND DATA 
ANALYSIS 

We selected a sample of clusters for which luminosities, cool- 
ing times and Chandra data were available. From there, 
we mostly limited ourselves to those that were cool cores 
(tcooi ^ 7 Gyr) and highly-luminous with average Lx ^ 
10^^ erg s~^ . The sample contains many of the luminous 
cool core clusters above a redshift of 0.08, but depends on 
what is available in the Chandra archive. Our aim was to 
select a large enough sample to show that our type of object 
is not rare, but not to select all clusters in the Chandra 
archive that met with the properties of our objects (i. e. 
clusters requiring strong feedback from their AGN but with 
no X-ray detectable nucleus). Our sample is therefore not 
complete. We suspect that there may be many other clus- 
ters that agree with the properties of our objects. 

For our selected sample, the first division of Table [1] out- 
lines the thirteen clusters that did not show any evidence of 
a point source in their Chandra images. Those in the second 
(Z3146, MS1455, R1347 and MS2137) show a hint of a point 
source. Finally, we also included the Ce ntaurus Cluster in or - 
der to compare our results with those of lTavlor et al.l (|2006|), 
as well as Cygnus A, since it is a perfect counter-examplqj 
with an obvious, but faint, point source. In summary, an 



^ Other examples of luminous clusters hosting a bright X-ray 
3oint source at their cent res include H 182 1+6 43 (|Russell et al.l 



120101 ) IRAS09104-f-4109 dlwasawa et al.ll200ll), MA CSJ1931.8- 
2634 dEhlert et al.ll20ldl V 3C295 (JAllen et al.ll200lD . 4C+55.16 
dlwasa w a et al.l 200lh , and even the Perseus Cluster (e. g. 
Fabian et al.ll2000l V 
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Table 1: Sample of Clusters - (1) Name; (2) redshift; (3) IR luminosity (Lir); (4) identification number of Chandra observational project; (5) total exposure time; 
(6) central cooling time in Giga-years; (7) unabsorbed bolometric X-ray luminosity within r = 200 kpc, with 2a upper and lower uncertainties; (8) cooling radius; 
(9) kinetic luminosity determined as the unabsorbed X-ray bolometric luminosity within r = 1.34rcooi, with 2a upper and lower uncertainties; (10) 2a upper limit of 
background corrected number of counts for the nucleus; (11) and (12) 2a upper limit count rate converted into an unabsorbed nucleus luminosity using pimms; (13) and 
(14) unabsorbed spectroscopic luminosities (rest-frame 2 — 10 keV) with an additional adsorption of the power law of A^h = (Column 13) and log(A^H) = 23 (Column 14) 
along with the 2a upper and lower uncertainties. When only upper limits were available, the 2a upper limit is shown. The dash symbol (-) indicates eit her no available 
data i n the literature for Column 3 or no fit to the data for Columns 13 and 14. Comments: ^- Galactic abs orption was f r ee to vary since the value from iKalberla et al.l 
(|2005|) did not provide a good fit. ^- The X-ray emission region for the nucleus was taken as the same one in lTavlor et al.l (|2006|). ^- Internal absorption of the power law 
was always allowed for Cygnus A where we used our spectroscopically derived value (log(A^H) = 23.33), but all luminosities shown are corrected for absorption. 
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absence of a detected point source is apparent in about half 
of the luminous clusters we examine. 

Table [T] includes the infrared (IR) luminosities of the 
BCG in our c lust ers. These luminosit ies were taken f rom 
iQuillen et all dJOOSi ) (8 - lOOO/xm), Egami et all (120061) f^ 
3 — TQ Ojtim) or using the 15 /xm luminosity of lEgami et al.l 
((20061) a nd the equat i on re lating the latter to the IR lumi- 
nosity of lElbaz et all (l2002l ) (8 - lOOO/xm). 

Table [1] shows the observation identification numbers 
and exposure times of the X-ray data. We took one of the 
deepest available non- grating data sets, except for Cygnus A 
where we decided to take OBSID 1707 which had a frame 
rate of only 0.4s to avoid pile- up ((Young et al.ll2002i ). For 
those that had data sets with similar exposure times, we 
selected the most recent one. The cooling times were taken 
either from Allen (2000), Bauer et al. (2005), Dunn & Fabian 
(2006), Rafferty et al. (2008) or Kirkpatr ick et al. ( 2009) . 
The cooling radius was taken from Dunn fc FabianI (|2006l . 
I2OO8I ) , or as the radius for whi ch the cooling time was equal 
to 3 Gyr in the online data of ICavagnolo et aU (|2009l ). For 
all of our objects, the radius where tcooi ~ 1 Gyr is just less 
than half the listed rcooi or 20 — 50 kpc. The region with a 
short (< 1 Gyr) cooling time therefore remains large. 

The X-ray Chandra data were processed, cleaned and 
calibrated using the latest version of the CIAO software 
(ciAOv4.2, CALD B4.2). Spectra were analysed using xspec 
(vl2.6.0d,f)(e. g. lArnaudlfToOel ). Fig.[T]shows images of four 
clusters in our sample. The top panels are images in the 
0.5 — 7 keV energy range, and the bottom panels are the 
high energy images (5 — 7 keV) where you would expect 
to have hints of a point source if a non-thermal component 
was present. Fig. [1] shows that the X-ray luminosity peaks 
toward the centre, but that there is no clear evidence of a 
luminous point-like AGN for R1504 and A1835. 

The cluster luminosities shown in Column 7 of Table 
[Uwere calculated by extracting a 0.5 — 7 keV spectrum of 
the inner 200 kpc, binned with a minimum of 30 counts per 
bin, and then fitt ing it with a single-temperature mekal 
(JMewe et al.l Il995l ) model including phabs absorption to 
account for Galactic absorption. The Galac tic absorptions 
throughout this paper were frozen at the iKalberla et al.l 
((2OO5I ) values, except for Centaurus and A0478 where the 
absorption was free to vary in order to provide a good fit. 
The redshifts were also kept frozen throughout this paper. 
Using the cflux model in xspec, we obtained a flux esti- 
mate for the cluster corrected for absorption, which was then 
converted into a luminosity using the luminosity distance. 
CFLUX provides a better estimate of the flux by exploring 
the space allowed by the free parameters. The kinetic lumi- 
nosities shown in Column 9 were calculated in the same way, 
but w ithin the inner region o f a radius (r) equal to 1.34 rcooi, 
where iDunn fc FabianI 



((2008[ ) found that on average bubbles 
can offset cooling within this radius. For both luminosities, 
the 2a upper and lower uncertainties are shown. All error un- 
certainties were found with the error command in xspec. 
The flux of the nucleus in each cluster was estimated 
using two methods, one photometric and another spectro- 
scopic. Thefirst was done using the web interface of pimm^j 
(|Mukailll993l ). For each cluster, the location of the nuclear 



http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 



X-ray emission in the 0.5 — 7 keV Chandra images was taken 
as the 4 by 4 pixels square at the brightest point in the cen- 
tre of the cluster. The background was taken as a square 
annulus centred on the same position, with an outer 8 by 8 
pixels square and an inner 6 by 6 pixels square. Since the X- 
ray counts are governed by Poisson statistics, the noise can 
be estimated as a = V^ in each region. The total count 
number of the nucleus X-ray emission was calculated by 
subtracting background emission scaled to the same num- 
ber of pixels. Using the error propagation equation, we then 
estimated the 2a noise level for the background corrected 
count number, and determined a 2a upper limit for the 
number of counts of the X-ray nucleus. This is shown in 
Column 10 of Table [1] Although background subtracted, 
this count number could still include some thermal contri- 
bution. The luminosities obtained are therefore upper limits 
for any non-thermal contribution. PIMMS was then used to 
convert the count rate into a flux and then luminosity in 
the 0.01 - 50 keV (Column 11) and 2-10 keV (Column 12) 
band. A power la w model with photon index of 1.9 was used 
(JGilli et al.ll2007l ). but our results are not sensitive to our 
value of the index, at least within ±0.2. For the Centaurus 
Cluster, we used the same location and region for the X-ra y 
nucleus emission and background as in lTavlor et al.l ((20061 ). 
in order to compare both studies. 

We investigated whether we could use the assumption 
that surface brightness scales with radius as a power law, 
since the surface brightness scales with the square of the 
density, and dens ity scales as a power law with the radius 
((Allen et al.ll2006l ). Averaging over each annulus, the count 
rate per pixel can be used as a measure of the surface bright- 
ness profile. The difference between any emission arising 
from the nucleus and the expected emission arising from a 
fitted power law to the surface brightness profile should give 
an estimate of the intrinsic nucleus emission. From this in- 
trinsic emission, pimms can then be used to estimate the flux 
of the nucleus. However, when fitting the surface brightness 
profiles to the inner ^ 10^^ with a power law, excluding the 
inner ^ 1", this predicted a higher number of counts than 
that found within the nucleus. Hence, our clusters have sur- 
face brightness profiles that flatten within the inner regions, 
and this extrapolation approach cannot be used to estimate 
the flux of the nucleus. 

PIMMS gives a first estimate of the luminosities based 
on a count rate, but it is useful to see if a non-thermal 
component is hidden within the spectrum of the nuclear X- 
ray emission. If there is a hidden component, this should 
provide a more precise estimate of the luminosity. The next 
two paragraphs describe in detail this method. 

First, we used the same location for the X-ray nucleus 
as that of the first method, but with a circular region of 
radius l^^ Second, another region was taken as an annulus 
centred on the same position, but with radii between 2" and 
3^^ Using this annulus as a background, and then fitting a 
power law model to the data did not provide a good fit. 
Instead, we used this annulus to estimate the properties of 
the cluster thermal component, which we then extrapolated 
to the inner circle of radius 1^^ (see next paragraph). The 
background was therefore taken as a region located far from 
cluster emission. Next, C-statistics were used to account for 
the low number of cou nts. However, we used the modified 
version of C-statistics (lArnaudI l2003l ) . which allows us to 



© 2010 RAS, MNRAS 000, ??^ 



Radiatively-inefficient cluster cores 5 



0.01 r 



10-^ 
2 



1.5 P 

1 




0.5 



W 10-3 r 



%w^ 



Energy (keV) 




Energy (keV) 

Figure 2. Example of one of our good spectral fits of the X- 
ray nucleus spectrum (A2204, top), and one of our worst fits 
(A1835, bottom). Each nucleus was fitted with a combined ther- 
mal (mekal, dotted line) and non-thermal (POW, dashed line) 
model. The combined model is shown with the black line. Galactic 
absorption was kept frozen at the value of Kalberla et al. (2005|), 
no internal absorption of the power law was allowed, and the 
background was taken as a region located far from the cluster. 
The data along with their error bars are shown with the black 
points, but were rebinned for plotting purposes until a minimum 
of 5cr and a maximum of 30 bins for A2204 and until a minimum 
of 4cr and a maximum of 20 bins for A1835. See Section 2 for 
more details on the fits. 



trapolated to the radius of the X-ray emission region (circle, 
r = 1"), using the re lation that T = ar^ where b ^ 0.3 
(|Voigt fc FabianI [2OO4I V The abundance is not expected to 
vary significantly from r = 1^^ to r — 2 — 3^\ Using this abun- 
dance and extrapolated temperature, along with a power- 
law index of 1.9, the only parameters allowed to vary in the 
model fitted to the nucleus emission (circle, r = 1'') were the 
normalization parameters of the thermal and non-thermal 
component. However, the normalization of the thermal com- 
ponent was not allowed to be less than the normalization ob- 
tained in the fit for the annulus (r = 2 — 3^^), scaled for the 
same pixel number. This is because the normalization pa- 
rameter is proportional to the square of the density which is 
expected to increase as the radius decreases. The additional 
internal absorption for the power-law was kept frozen at ei- 
ther A^H = (Table 1, Column 13) or \og{Nu) = 23 (Table 
1, Column 14). See also Fig. [2] for some examples. Finally, 
we extracted the unabsorbed flux of the non-thermal com- 
ponent in the 2 — 10 keV band using the CFLUX model. The 
2a upper and lower uncertainties are also shown in Table 1, 
but some clusters only yielded upper limits, as shown by the 
< symbol followed by the 2a upper limit. 

We stress that we were able to obtain a rough estimate 
of the nucleus luminosities based on the spectral analysis, 
but required most parameters to be frozen. It was not clear 
whether the addition of the power law model did actually 
improve the fit. Since our analysis is based on C-statistics, 
it is not possible to provide a goodness of fit like that of x^ 
statistics. The lack of any obvious non-thermal contribution 
to the nucleus spectrum of our objects supports the idea 
that they are radiatively-inefficient. 

X-ray studies of the innermost regions of nearby cool 
core clusters typically show several cooler components. We 
strongly suspect that the apparent positive detections of 
non-thermal components in the present spectral analysis are 
just manifestations of these cooler thermal components. We 
therefore proceed by treating all such nucleus detections as 
upper limits. 



use a background file by loading it in xspec with the back 
command. We then fitted the spectrum of the nucleus (circle, 
r = 1^^) in the 0.5-7 keV band with a mekal + power law 
model in xspec, including Galactic absorption (phabs). Ad- 
ditional internal absorption (zphabs) of the power-law was 
allowed. However, the power law component made the model 
largely unconstrained. Even for the clusters which had long 
exposure times such as A2204, it was difficult to obtain any 
constraints on the power law parameters. In order to provide 
a good fit, it was therefore necessary to freeze most parame- 
ters. The next paragraph outlines the parameters that were 
kept frozen, and explains how the fit was accomplished. 

First, the Galactic absorption was kept frozen. Next, 
and as mentioned earlier, the spectrum of the second region 
(annulus, r = 2 — 3^^) was used to estimate the properties 
of the cluster thermal component. We fitted an absorbed 
(Galactic) single-temperature mekal model to the annulus 
where only the temperature, abundance and normalization 
were free to vary. The extracted temperature was then ex- 



In the case of Cygnus A, the spectrum clearly showed 
a non-thermal component which required a specific internal 
absorption. The spectroscopic method was first applied to 
estimate the luminosity of the nucleus, in the same way as 
for the other clusters but where the internal absorption of 
the power law was free to vary. The fitted model gave a 
value for this absorption of 21.2t\\l x 10^^ cm"^ (log(A^H) = 
23.33). Column 14 shows the derived luminosity corrected 
for absorption. Using this additional absorption, the nucleus 
luminosity was then estimated using pimms in the 0.01 — 
50 keV and 2 — 10 keV bands, see respectively in Column 
11 and 12. The unabsorbed cluster and kinetic luminosities 
shown in Column 7 and 9 include the contribution of the 
non-thermal component of the nucleus. Here, the internal 
absorption was also kept frozen to the value derived for the 
nucleus (log(ArH) = 23.33). 

Projection effects cannot be corrected for accurately 
enough near the central regions of clusters, and for this rea- 
son we did not consider them throughout this paper. 
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3 RESULTS 



The r esults are essentially shown in Table \T\ iTavlor et alJ 
(|2006l ) found a 3a upper limit of 1.2 x 10"^° erg s"^ for 
the nucleus of Centaurus. They used the annulus surround- 
ing the nucleus as a background and by combining vari- 
ous data sets (200 ks) were able to constrain more pa- 
rameters including the internal absorption of the power 
law (A^H '^ 3 X 10^^ cm~^). We compare this with our 
PIMMS method which also uses the surrounding annulus 
as a background. Including the same absorption and our 
2(7 upper limit count rate, we derive an unabsorbed upper 
limit of 1.56 X 10^° erg s~^ for the nucleus, consistent with 
iTavlor et all (| 20061 ). 

For Cygnus A, we find an internal absorption of 



21.2; 



-1.5 

-1.4 



^22 



X 10 cm and a 2 — 10 keV unabsorbed luminos- 



ity for the nucleus of 2.11q'2 x 10^^ erg s ~^ , using a power- 
law index of 1.9. For the same object, lUeno et al.l (| 19941 ) 
found an internal absorption of 37. Sl^^ x 10^^ cm~^ and a 
2—10 keV unabsorbed luminosity of ^ 5.5 x 10^^ erg s~^ 
(using our cosmology). Their best-fit power-law index was 
1-98-0.20 5 ^^^ they estimate at the 90 per cent confidence 
level that their intrinsic luminosities are uncertain by a fac- 
tor of 2. Within this uncertainty, our results are consistent 
with theirs. Our r esults are also consistent with those of 
I Young et al.l ((2003). These authors used the same Chandra 
data as our study, and found a photon index of ^ 1.5, an 
internal absorption of 20^2 x 10^^ cm~^ and a 2 - 10 keV 
unabsorbed luminosity of '-^ 1.9 x 10^^ erg s~^ (using our 
cosmology) . 

In order to illustrate the power needed to offset cool- 
ing and the radiative power o f the AGN, as w e ll as com- 
pare our results with those of iMerloni fc Heinj (|2007l V we 
plot the kinetic power versus the nuclear luminosity in the 
2-lO keV band (Fig.[3l). Here the results of lMerloni fc Heind 
([20071) are shown with black filled circles. Using the par- 
tial Kendall's r correlation test, they found that the kinetic 
power was correlated with the nuclear X-ray luminosity to 
more than a 3a level. In order to guide the eye, we plot in 
Fig. [3] the linear regression to these data in a log-log plot 
with a black line, along as the 2a upper and lower limits 
(dotted lines). The linear regression consisted of minimising 
the x^ 5 cind the error bar estimates of Lkin were taken as the 
mean value between the upper and lower error bar measure- 
ments. iMerloni fc Heinj (|2007l Vs data included two objects 
which only had upper limits of Lnuc- Our linear regression 
does not include error estimates of Lnuc- We therefore omit- 
ted these two objects for the regression. The non- filled sym- 
bols in Fig. [3] illustrate our results. Here, the kinetic lumi- 
nosities and upper limit of the nucleus luminosity are taken 
from Table [U the later using the upper limit values derived 
with PIMMS. For Cygnus A we used the spectroscopic derived 
value of Lnuc, since it is clearly detected. The clusters that 
have square symbols also have detailed estimates of their jet 
power, obtained from observed bubbles. For R1532, A1835, 
A2204, A0478, P0745, Z2701 and Z3146, this estimate is 
shown with the dark blue filled squa r e and w as taken in order 
of preferenc e from [Rafferty et al.l (|2006l l: iDunn fc FabianI 
([2008 ); Sanders et al.l (|2009|). Cygnus A is shown in green, 
as well as its jet power estimate in filled green taken from 
IMerloni Sz Heinj (|2007l l. Centaur us is shown in red an d its 
jet power in filled red taken from iRaffertv et al.l ((20061). 
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Figure 3. Kinetic bolometric luminosity versus nuclear luminos- 
ity in the 2-10 keV band. The AGN of Merloni & Heinz (2007) 
are shown with black filled circles. Our results are shown with 
the non- filled symbols, where our value for Lcool(''^ ^ l-34rcool) 
is used as the kinetic luminosity and the upper limit PIMMS esti- 
mate of the 2 — 10 keV nucleus luminosity is used for Lnuc- For 
Cygnus A, since its power law component is clearly detected, we 
used our derived spectroscopic value of Lnuc- The green symbols 
are for Cygnus A and the red symbols are for Centaurus Cluster. 
The clusters that have estimates of their kinetic luminosity from 
measurements of the energy output from bubbles/jets are shown 
in squares and the estimate is shown with the filled square. The 
black line shows the linear fit to the data, and the dotted lines 
are the 2a limits. 



This figure shows that our sample of clu s ters m ostly lie 
well above the relation of Merloni fc Heinj (|2007l ) 's data, 
but that Cygnus A and Centaurus fall well inside the con- 
fidence levels. All of the clusters in our sample that do not 
have an obvious X-ray point source have jet powers above 
lO^^"^ erg s~^ . Hence, they seem to require a jet power from 
their AGN that is very high. On the other hand, if our ob- 
jects are simply shifted to the left because they are highly 
obscured, they would require a large intrinsic absorption of 
about log(ArH) = 23.6. 



4 DISCUSSION 

Our results reveal a significant population of objects re- 
quiring high kinetic input from an AGN to offset cool- 
ing and/or high jet power, yet are without a detected X- 
ray nucleus. These objects appear to be radiatively ineffi- 
cient with, on average, Lkin/ Lnuc exceeding 200. Our clus- 
ters have radio luminosities at 1.4 GHz of about 10^° — 
10^^ erg s~^ , but require very strong jet powers to offset 
cooli ng (> 10^ ^ '^ er g s~^ ). When comparing these results 
to Birzan et al.l (|2008 ) 's relation between jet (cavity) power 
and Lr at 1.4 GHz, we find that our clusters lie above this 
relation, implying once more than our objects require ex- 
treme kinetic feedback. 
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4.1 Stellar mass black holes 

The upper-hmits to the X-ray luminosity of the nucleus in 
our sample are at least 10 times lower than is straight for- 
wardly predicted from the pro perties of the detected lower 
jet power objects discussed bv lMerloni fc Heinj (|2007|) (see 
Fig. [3|). Stellar mass black hole systems may have similarly 
high ratios of kinetic to luminous power, such as Cygnus X- 
1 which is thought to spend most of its time (90 per cent) 
in a low state, w ith a hard spectrum and collimated jet. 
iGallo et al.l (|2005l ) found that the jet in Cygnus X-1 appears 
to be inflating a large ring-like structure, and that the kinetic 
power required can be as high as the bolometric X-ray lumi- 
nosity of the system. These authors then suggested that, in 
general, stellar mass black holes dissipate most of their ac- 
creted power into radiatively-inefficient relativistic outflows, 
and not in the form of X-ray emitting i n flow. A similar con- 
clusion was reached by iKording et al.l (| 20061 ) but for BHs 
also found in AGN, i.e. that a BH in a hard state looses 
most of its energy in jet power. iKording et al.l (|2008|) then 
found that all low- luminosity AGN seem to launch strong 
jets. 

We do not know the black hole (BH) mass for most 
of our objects but we expect that they lie between 10^ and 
10^° M0, based on other Brightest Cluster Galaxies (BCGs). 
Given that these objects are in an exceptional environment 
and are active continuously, it is not clear that the standard 
Mbh — ct or Mbh — Mk relations are relevant for them. The 
few BCGs with a rehable Mbh (e. g. M87, Macchetto et al. 
1997, as weU as A3565, A1836 and A2052, Dalla Bonta et 
al. 2009) indicate the above mass range. The kinetic power 
the n corresponds to 1 to .1 per cent of the Eddington limit. 

IChurazov et al.l (|2005h have argued, by analogy with the 
interaction of X-ray binaries, that black holes operating at 
powers close to the Eddington limit (LEdd) are radiatively 
efficient while at low power (less than 1 per cent LEdd), 
their efficiency drops steadily to low values. The bulk of the 
power is increasingly taken up by outflows at low Eddington 
rates. This explains the behaviour of the low and intermedi- 
ate luminosity galaxies, groups and clusters but is difficult to 
accomm odate for the highest lu minosity objects found here 
(see also iMcNamara et al.ll2010l ). The power emerging from 
our objects means that they are operating above 10~^LEdd 
if Mbh ^ lO^M©. The radiative power should then be at 
least comparable to the kinetic power. However, if the cen- 
tral black holes significantly exceed lO^^M©, i. e. are ultra- 
massive, then the radiative power could fall well below the 
kinetic power, which is conceivable for our powerful objects 
but remains untestable at present. 



4.2 Relativistic and geometric effects 

Among the possible explanations as to why our objects have 
such low observed X-ray luminosities, we first mention that 
they might simply be Doppler-suppressed. If the X-rays are 
mostly produced in the jet, and if the jet axis of the AGN 
lies close to the plane of the sky, as well as if the jet has a 
high Lorentz factor, the core luminosity would be suppressed 
by a significant factor such that Lnuc = D^ x Lnuc,o- Here, 
Lnuc is the observed value of the core luminosity, Lnuc,o is 
the intrinsic luminosity and D is is the Doppler factor given 
by D = l/(r(l - /3 X cos 6)) (F is the Lorentz factor and 



is the angle between the line of sight and the jet axis). 
However, in order for o ur objects to sit on the relation of 
iMerloni fc Heinj (|2007|), their jet axis would not only have 
to lie close to the plane of the sky, but the jet velocity would 
also have to be such that /3 ^ 0.9. It would be unlikely that 
all of our objects have jets with a preferred geometry (axis 
aligned with the plane of the sky) , and since our sample is 
not statistically complete, it is possible that there are other 
objects that meet the properties of those in our sample. 



4.3 Advection dominated accretion floAvs or 
magnetically-dominated black holes 

Second, we mention that they might be standard advec- 
tion dominated accretion flows (ADAFs). Although ADAFs 
models can create strong winds and relativistic jets, it is 
not certain that they would be able to create the extreme 
relativistic jet powers required to inflate the bubbles in 
our objects, since ADAFs los e much matter in their winds 
(JNaravan fc McClintoc3l2008l ). 

We also consider whether the y are powered by 
magn etically-dominated accretion (e. g. iBlandford fc Pavnd 
Il982i ). Our objects have so far shown no evidence for non- 
relativist ic winds. They may instead be simply magnetically- 
dominated systems (more than 90 per cent) harbouring 
powerful jetted outflow s. iKuncic fc Bicknelll (|2007|) and 
iMerloni fc FabianI (|2002h pointed out that energy and an- 
gular momentum could be removed from an accretion disc 
in the form of powerful Poynting-dominated outflows if 
strong magnetic fields were present, while being radiatively- 
inefficient. The mechanism responsible for creating jets in 
magnetically-dominated accretion discs still remains poorly 
understood, which makes it difficult to give any definitive 
conclusion. 



4.4 Highly-absorbed AGN 

On the other hand. X-ray absorption can render the nucleus 
undetectable in 0.5 — 7 keV X -rays. The power m ust then 
emerge at longer wavelengths. lEgami et al.l (|200Q) studied 
Spitzer data of a sample of eleven brightest galaxies in X- 
ray luminous clusters. This sample included Z3146, A1835, 
MS1455, A2261, MS2137 and A2390. Essentially, they found 
that these clusters, which had the shortest cooling times in 
their sample, had the largest BCG IR luminosities. This 
seems to suggest that highly-luminous strong cool core clus- 
ters have BCGs with larger IR luminosities. These authors 
also found that A1835 and Z3146 could be classified as 
LIRGs. Large IR luminosities could be a sign of strong ab- 
sorption, which would render the X-ray nucleus almost in- 
visible^ 

ICandhi et al.l (|2009h reported that there was a 1:1 cor- 
relation between the intrinsic 2—10 keV X-ray luminosity 
of an AGN ( Lnuc) and its mid-IR luminosity at ^ 12/im, 
impl ying that both quantif ies are intrinsically related (see 
also iHardcastle et al.l |2009|, who looked at a more com- 
plete sample of radio sources). If we assume that our objects 
must follow this relation (although they lie at the centre of 
some of the most extreme cool core clusters), then those 
that have large mid-IR luminosities (i.e. R1532, A1835, and 
Z3146) must have large intrinsic nuclear X-ray luminosities 
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on the order of 10^^ erg s ^ . This would mean that they 
are highly obscured, since they are almost invisible in our 
Chandra data with Lnuc(2 - 10) keV - 10^2-43 ^^^ ^-i _ 
The intrinsic nuclear X-ray luminosities, as predicted from 
the X-ray /mid-IR relation, w ould also make them consistent 
with lMerloni fc Heind (|2007l Vs relation (see Fig.E)). For the 
remaining objects in our sample, if they are required to fol- 
low the X-ray /mid-IR relation, then this would imply that 
they are not highly obscured, since the predicted Lnuc would 
be consistent with those that we find (i. e. ^ 10^^ erg s~^ ). 
It remains that with the advent of ALMA (Atacama Large 
Millimeter/submillimeter Array), we will be able to resolve 
in detail the dusty structures within the core of clusters 
through sub-mm observations, which will allow us to con- 
firm if the large IR luminosities of some of our BCGs are 
of nuclear origin. We must stress, however, that most of the 
far IR emission in our objects is likely due to young stars. 
There is much evidence from excess blue ligh t, that there 
is on going star formation in these objects (see lO'Dea et al.l 
[2OO8I). Nevertheless, there is still room for a significant AGN 
contribution. 

If our objects have different dust properties than 
those that f o how the X-ray/mid-IR relation found by 
iGandhi et al.l (|2009|), they could still be heavily absorbed. 
In Column 14 of Table [1] we have shown that even with 
a log(A^H) = 23, resembling the absorption of Cygnus A, 
the X-ray nucleus luminosity is at most 3 times higher than 
without internal absorption. If our objects are very heav- 
ily absorbed (log(A^H) = 24), this would make their non- 
thermal component almost invisible in the 0.5 — 7 keV X- 
ray energy band. However, our objects would still represent 
a different population from Cygnus A which has an obvi- 
ous point-like core accounting for more than 50 per cent of 
Lx even with strong absorption (log(A^H) = 23.33). With 
today's high energy X-ray telescopes, it remains difficult to 
properly image the high energy band, since the background, 
sensitivity and spatial resolution remain problematic. How- 
ever, in the near future, satellites such as NuSTAR (nuclear 
spectroscopic telescope array, 6 — 79 keV) and Astro — H 
(0.3 — 80 keV) could provide the sensitivity and resolution 
needed to analyse the power-law component of our objects 
at high X-ray energy. 



4.5 AGN duty cycles 

Fifth, we consider strong variability. The duty cycle of an 
AGN, which measures the probability that a BH is in an 
active phase is less than 10 per cent (e. g. IShankar et al.l 
[2OO9I). For our objects, this requires that the power be above 
10^^ erg s~^ when the AGN is swi tched on. If the duty cy - 
cle is less than 1 per cent (e. g. ICiotti fc Ostrikerl 120071 ) . 
the power must exceed 10^^ erg s~^ , which is greater than 
the most luminous jets known and would surely disrupt the 
core. Such observed high power jets are not associated with 
BCGs, but rather blazars. T he only quasar in a BCG be- 
low z = 0.5 is H1821+643 (iRussell et al.ll2010l ). In other 
words, if this explanation is correct, then contrary to ob- 
servations there should be a population of powerful quasars 
and blazars at low redshift s associated with BCGs. How- 
ever, iDunn fc FabianI (|2008r ) found that out of a sample of 
34 clusters, only 6 had extreme outbursts capable of offset- 
ting cooling to more than twice the cooling radius. Of these. 



Z3146 and Z2701 were included. This implies that the duty 
cycle of extreme outbursts is about 18 per cent, which would 
require that the power be about 5 x 10^^ erg s~^ when the 
AGN in our sample are in this mode. 

On the other hand, it is also possible that the proper- 
ties of our objects are not related to the variability of the 
jet itself, but rather the variability of how many X-rays are 
produced by the core, i.e. the radiative efficiency of the jets. 
Our objects may have jets that are continuously pumping 
10^^ erg s~^ into the medium (which would provide the nec- 
essary amount of energy needed to heat the ICM) , but are 
in a very radiatively inefficient mode. A scenario here could 
be to assume that the radiate efficiency of a jet depends 
on its environment, i.e. depends on just how much interac- 
tion directly occurs between the jet and surrounding gas. 
Extremely powerful steady jets may create a channel which 
results in little such interaction for the bulk of the jet. 



4.6 Spin-poAvered black holes 

Finally, our black holes could also be powered by spin, not 
accretion (e. g. iBegelman et al.l ll984l V iMcNamara et al.l 
(|2009l ) studied the case of MS0745.6+7421 harbouring one 
of the most powerful outbursts ('^ 10^^ erg) but which has 
no X-ray point-like core in the Chandra image. They sug- 
gest that the outburst could be powered by angular mo- 
mentum originat i ng fr om a rapidly rotating BH (see also 
iMcNamara et al.l I2OIOI V Strong jets could be powered by 
rotating BHs (see e. g. the BZ mechanism of Blandford 
& Znajek 1977, but see also Fender, Gallo & Russell 2010 
for a discussion on the observed lack of any correlation be- 
tween jettedness and spin in galactic microquasars) . To op- 
erate a spin-powered feedback loop and maintain the bal- 
ance between cooling and heating, the energy must be con- 
tinuously tapped. A high spin parameter is required, which 
according to our significant population of objects would re- 
quire that highly rotating BHs are not rare. In order to test 
this theory indirectly, one could look at the relation found 
by Allen et al. (2006). In their sample of nearby ellipticals, 
which includes galaxies with an obvious point-like X-ray nu- 
cleus, the Bondi accretion rate correlates with jet power. If 
the spin explanation is correct, then this correlation would 
not hold for all objects. 



5 SUMMARY 

We have identified a population of objects which require 
powerful jets to be present but have no X-ray detectable 
nucleus. We have also identified a range of possible expla- 
nations, each of which carries significant implication for the 
origin and operation of jets. The black holes may be ul- 
tramassive (Mbh ^ lO^M©), or have very high spin, or be 
highly obscured. They may also be mostly off, yet unobserv- 
able when they are switched on, or have highly radiatively 
inefficient jets. 
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